Continuous phase transition in a spin-glass model without time-reversal symmetry 
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We investigate the phase transition in a strongly disordered short-range three-spin interaction 
£vq \ model characterized by the absence of time reversal symmetry in the Hamiltonian. In the mean- 

field limit the model is well described by the Adam-Gibbs-DiMarzio scenario for the glass transition; 
QQ ' however in the short-range case this picture turns out to be modified. The model presents a finite 

^\ , temperature continuous phase transition characterized by a divergent spin-glass susceptibility and 

a negative specific heat exponent. We expect the nature of the transition in this 3-spin model to be 
the same as the transition in the Edwards- Anderson model in a magnetic field, with the advantage 
that the strong crossover effects present in the latter case are absent. 

T3 

£h , I. INTRODUCTION 

o 
o 

Nowadays there is large activity of research in the problem of the glass transition from the perspective af spin-glass 
theoryBa. This interest originates from old observations by T. Kirkpatrick, D. Thirumalai and P. Wolynesta who found 
a striking similarity between the dynamical equations of some mean-field spin glass models and the mode-coupling 
equations for glasses. The mode-coupling equations are characterized by the presence of a dynamical singularity at a 
temperature Td below which spin-spin correlation functions do not decay to zero in the infinite time limit signalling 
the breaking of ergodicityu. Above but close to Td the correlation functions display a plateau which separates two 
different relaxational regimes (the alpha and the beta processes) . The family of models which show this behavior are 
those with one-step of replica symmetry breaking (models with an infinite number of breaking steps describe better 
the spin-glass behavior found in strongly disordered magnets). These models are characterized by two singularities 
or transition temperatures. One transition is purely dynamical and corresponds to the mode-coupling transition Td 
previously described. The other transition at T c < Td is thermodynamic and corresponds to a temperature below 
which replica symmetry breaks and the configurational entropy (also called complexity) vanishes. The transition at T c 
has features of both first and second order transitions like a discontinuity in the Edwards- Anderson order parameter 
and a finite jump in the specific heat. 

It is widely accepted that the origin of the dynamical transition at T<j (where the relaxation time diverges and 
ergodicity breaks) originates from the presence of an exponentially large number of states (exponentially large with 
the system size) which trap the system for exponentially large times forbidding the system to reach the equilibrium 
Gibbs measure. But it is also clear (and this was also recognized as a strong limitation in the original mode coupling- 
theory) that equilibrium below Td should be restored in finite-dimensional systems where activated or nucleation 
process (i.e. jumps over finite free energy barriers) take place in a finite time. What is the final theory which correctly 
describes the nucleation processes taking place in glasses is one of the major open problems. The behavior of T c in 
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the presence of short-range interactions is less clear. According to the mean-field picture the transition at T c (where 
replica symmetry breaks) could well survive in finite dimensions. At T c the configurational entropy would still vanish. 
This scenario is nothing else (but now rephrased in the spin-glass language) than the Adam-Gibbs-DiMarzio scenario 
(hereafter referred as AGM) for the ideal glass transitiontm. 

On top of this connection between the spin-glass theory and mode-coupling theory, a very interesting connection has 
also been established between the statics and dynamics of glassy systems in the off-equilibrium regime. In this case, 
the equilibrium order parameter for spin-glasses (the so called P(q) function) is intimately related to the fluctuation- 
dissipation ration. Thi&Jink between statics and dynamics, originally suggested by the analysis of exactly solvable 
mean-field spin glassesQa, is actually supported by extensive numerical simulations!! and general arguments based on 
the assumption of linear response theory applied to short-range spin-glass modelaid. 

All these previous studies cannot predict how short-range corrections change the mean-field behavior. And in 
particular, it is unclear how the AGM scenario typical of a first order spin-glass transitions is modified in finite 
dimensions. The answer may crucially depends on the presence of quenched disorder in the system and the mean-field 
critical behavior can be significantly altered in finite dimensions. Here, we will see that the mean- field scenario a la 
AGM does not survive in finite dimensions for a certain class of models. |— ■ 

A preliminary account of some parts of this work has recently appeared in a different contexttil where a Binder-like 
parameter was proposed to study replica rSymmetry breaking transitions. Also, different cases of the present model 
have been already studied in referencestHilil. So our work complements these results although the case we study here 
lacks time reversal symmetry sharing some features of the Ising spin glass in a magnetic field. 

The paper is divided as follows. In the section II we define the model and the numerical algorithm is explained 
in section III. In section IV we present the equilibrium results obtained by simulating small systems. In section V 
we present results for the order parameter and its cumulants. In section VI we give results for a new parameter 
(to be defined latter) which unambiguously shows the existence of a phase transition. Section VII is devoted to a 
discussion of the dynamical properties. Finally, in section VIII, we summarize the main results and discuss some 
peculiar features of the present model. 



II. THE SHORT-RANGED P-SPIN MODEL 



The, present model is a short-range generalization of the multi-spin interaction mean-field Ising spin glass defined 

bvLM 

H = - ^ Jiri2-i P < J ii (J i2--- (J i P > (!) 
(ii,i 2 ,..,i p ) 

where the spins (Tj, I < i < N (N is the size of the system) can take the values ±1 and the Ji 1 i 2 ..i p are quenched 
random variables with zero mean and variance p\/(2N p ~ 1 ). In ([!]) all possible multiplets of p spins interact through 
the random couplings Ji 1 i 2 ..i p ■ Consequently there is no spatial dimensionality and the model retains its full mean-field 
character. In order to go beyond mean-field theory we need to suitably modify the model introducing short-ranged 
interactions in a finite-dimensional lattice. A possible way to modify eq.(|l|) is to consider only links which couple 
nearest-neighbors in a finite-dimensional regular lattice. For instance, we could locate the spins in the vertices of 
the lattice and consider only a certain set of triangles (j> — 3), axectam set of squares (p = 4) and for a general 
p, only a certain set of plaquettes containing p spins-.aLthe latticeEZHlfl. In- the presence of quenched disorder, such 
constructions have been considered in several casesEJrEJ. In particular, ir£3 a simple cubic lattice with p = 4 was 
studied. Although freezing behavior was observed, no evidence for a finite temperature transition was found in three 
dimensions. Unfortunately, this type of models with binary exchange couplings has a large ground state degeneracy 
causing strong crossover effects (due to the presence of a non trivial zero-temperature fixed point) at low temperatures. 
Furthermore, in this type of models the lower critical dimension seems to increase with p. This implies that large 
dimensions need to be studied in order to find a phase transition. More work is certainly necessary to identify whether 
the AGM scenario is valid in finite dimensions for this type of models. 

The model we are going to study here is an alternative way to include short-range corrections. We locate M different 
Ising spins (J y i2 '" , ' D in each site of a regular cubic lattice. In this notation j enumerates the different spins (it ranges 
from 1 to M) in a given site with coordinates 12, in) where 1 < i\, Z2, «d < L and L is the lattice size of the 
cubic lattice and D is the dimensionality. The volume of the system is therefore given by V — L D . The Hamiltonian 
is defined as follows : 

D 

"H = ^ ^ T~tlink ) (2) 
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where TLu n k is the Hamiltonian corresponding to the link defined by the site («i, «2, and the direction fi, 

1 < fj, < D. In our notations, a link is a pair (P, /i) which couples the point P = (ii, *2, *_d) to the nearest neighbor 
site in the \i direction P + /i = (ii, .., i M + 1, i£>). Note that in eq.(||) each link is counted only once. For each link 
we sum all the possible groups of p spins out of the 2M spins located at nearest neighboring sites of the lattice (with 
p < 2AI). The final expression for TLu n k is given by 

p 

nj \ "* \ "* \ "* T (a 1 ...a k ,b 1 ...b p _ k ) p p p P+p, P+H P+» fo\ 

fe=l l< Qi <A/ l<fci<M 

The couplings j(" 1 -"°*' 6i "- 6 j , -»>) are random variables (which take the values ±1) uncorrelated for different links (P, /i) 
and sets ofp spins [a\...ak, bi...b p -k). Other versions of the model (for instance, ^f-*) = j(o l ...ofc,6i...b,_ k ) ) 

i.e. translational invariant disorder) are also possible and they could have different properties. 

As we observed in the introduction, the preserd-model has received considerable attention quite recently. A prelim- 
inary short account of our work was presented inlid. In three-dimensions, an exhaustive numerical study of the statics 
and the dynamics has been done in tha.case M = A,p = 4t3 while the case M = 3,p — 4 has been studied ir£3. So 
the main results on this model (exceptEJ) were obtained for the p = 4 model. The study pf the Gaussian propagators 
around the mean-field limit M — ► oo as well as the 1/M expansion were considered inEJ. The case we study here 
has the crucial property that the Hamiltonian (^) does not have the symmetry under time reversal (i.e. the global 
symmetry Oi — > — <7j, Vi). 

Here we will present exhaustive results for the case M = 2, p = 3 in D = 4. We have chosen these parameters for 
the following reasons: 

• p = 3: This is the simplest case which lacks time reversal symmetry. We expect in this case more clear results 
about the existence of replica symmetry breaking transitions in short-range systems. 

• M = 2: This is the simplest and non-trivial case. Larger values of M require always more computational effort. 

• D = 4: To be sure that we find a finite-temperature phase transition we have studied a large dimensionality 
compatible with a reasonable computational effort. 

For p = 3, M = 2 the general Hamiltonian eq.(||) reads 

V D 

IT. = - 2^ 2^ J (ll,10) Cr l Cr 2< J i + J(ll,01) Cr l Cr 2 cr 2 

l—l fJl— 1 

' J (10,11) CT 1 CT 1 "2 ' J (01,H) 2 1 2 ' ' ' 

where the (e^; /i = 1, .., D) denote the different unit vectors in a D dimensional lattice. The J's are binary uncorrelated 
random variables and we will consider periodic boundary conditions. 

Note that the model (||) has three spin interactions. One would simply expect the transition to belong to the class 
of (j) 3 theories. As we will see in the following sections, there is indeed a phase transition occurring in the present 
model. The absence of time-reversal symmetry in this model has crucial implications on the type of phase transition. 
We anticipate that the transition is related to the breaking of ergodicity at low temperatures, a consequence of the 
breaking of replica symmetry, the crucial symmetry to describe strongly disordered systems. We will try to clarify 
and give evidences on this point in forthcoming sections. 



III. THE NUMERICAL ALGORITHM 

We have studied the model in four dimensions using the parallel tempering methodil'El This is a good numerical 
method to equilibrate disordered systems at low temperatures. Contrarily to the simulated tempering method, in 
this algorithm it is not necessary to determine the free energy at different temperatures to reach equiprobability 
in the occupancies of these temperatures. Although the parallel tempering is a very efficient method to surmount 
energy barriers it is not clear how good the performance of the algorithm is in presence of entropy barriers (i.e. when 
relaxation to equilibrium takes place along narrow channels or gutters in phase space). 

The implemeatation of this algorithm is quite easy. It has been widely explained in the literature (for instance, 
see the reviewsc3'c3c3) and we will limit ourselves to sketch the main steps of the algorithm. We consider a set of 
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Nt copies or replicas of the same system which stay at different temperatures (Tj; i = 1,..,Nt). Each copy or 
replica is then specified by a pair (C, i) where C denotes the microscopic configuration (i.e. the values of all the spins) 
and the temperature of the copy i. We can then construct a Markov process in the space of configurations plus 
temperatures which satisfies ergodicity and detailed balance by allowing the following moves: 

1) Change configuration at fixed temperature: (C,i) — > (C',i) with probability 

P((C, i) -» (C, £)) = Mm(l, exp(-/3(«(C) - W(C))) . (5) 

2) Exchange configurations of two systems at temperatures (3i,(3j: {(C, i), (C ,j)} — > {(C, i), (C, j)} with probability 

p({(C,i),(C',i)}^{(C',i),(C,i)}) = M*n(l,exp(-(ft-&)(«(C)-«(C')) . (6) 

The first type of move is the usual change of configuration at fixed temperature in the Monte Carlo method. The 
second move avoids the system to get trapped in deep metastable minima. With this algorithm configurations which 
are far from each other can be reached by allowing a single copy of the system to extract energy from other copies 
through a coupling mechanism step (2), eq. (^) induced by the dynamics itself. Plainly speaking the rest of the copies 
or replicas play the role of an external bath for a given copy. If a given copy remains trapped in a deep minima 
of the free energy it can escape by extracting energy from the rest of the copies. The full Markov chain reaches 
thermal equilibrium when all temperatures are equally occupied and for each temperature the conditioned probability 
distribution P(C\i) is a Boltzmann distribution at temperature l//3j. In this way one simulates the model at different 
temperatures in the same run while being always in thermal equilibrium at different temperatures. 

The time needed for the Markov process to reach the stationary distribution (i.e. the thcrmalization time) depends 
mainly on the choice of the set of temperatures and also on the ratio between the number of moves of the first (§) 
and second type @. To solve the first problem, we have chosen a set of temperatures equally spaced in j3 = 1/T 
in such a way that moves of the second type (^J) do not occur with a too small probability at low temperatures. 
On the other hand, to uncorrelate the configurations as much as possible (in order to explore maximally distant 
configurations) it is convenient that the copies reach high enough temperatures. This is accomplished by enlarging 
the set of temperatures in the simulation up to twice the value of the critical temperature. For the present model 
we simulated 25 temperatures ranging from T = 2.0 up to T — 5.0 (as we will see this corresponds to a window of 
temperatures covering the region 0.75T C — 2T C ). As previously remarked, a good thcrmalization can be achieved by 
choosing an appropriate number for the ratio between the number of tempering moves of the first and second type. 
If the number of moves of the second type is too small then the system is not able to efficiently decorrelate and 
only a small number of excursions into the high-temperature region are performed. In the other extreme, the system 
decorrelates too fast and does not efficiently sample the landscape at a given temperature. We have tried several 
schemes intermediate between these two extremal cases and we have found that one move of the second type for each 
ten moves of the first type is a good compromise which efficiently thermalizes reasonable sizes. 

We did extensive simulations for L = 3, 4, 5, 6 with 2 15 = 32768 MCS (L = 3, 4) and 2 18 = 262144 [L = 5, 6) MCS 
at each temperature. For sizes L = 3, 4, 5, 6 we studied 1000, 600, 300, 100 samples respectively. A preliminary study 
for a set of 5 samples showed that this number of steps was enough to reach thermal equilibrium for the selected sizes 
for the range of temperatures studied. This preliminary study turned out to be crucial to determine the smallest 
range of temperatures which allows the system to equilibrate. If the range of temperatures selected extends down to 
too low temperatures then thcrmalization is hardly achieved. This means that the full Markov process associated to 
the parallel tempering algorithm does not reach the stationary Boltzmann solution. Then, there is no guarantee that 
thermalization is achieved neither at low nor at high temperatures. 



IV. THERMODYNAMIC OBSERVABLES 



A preliminary research of the evidences for a phase transition includes the study of the temperature behavior of 
extensive quantities such as the internal energy or the specific heat. In figures 1 and 2 we show the internal energy 
and the specific heat averaged over the samples for different sizes as a function of temperature. The energy E and 
the specific heat C were computed using the expression, 

E = Wt (7) 

c = ^{W)-W?)- (8) 
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2 3 4 5 
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FIG. 1. Energy versus temperature (L = 3,4,5,6 correspond to long dashed, short dashed, dot and solid lines). Error bars 
are shown for L = 6. 
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FIG. 2. Specific heat versus temperature (L — 3,4,5,6 correspond to long dashed, short dashed, dot and solid lines). Error 
bars are shown for L — 6. 

In what follows < .. > stands for Gibbs average and (..) for disorder average. We note the absence of any jump in 
the internal energy as well as divergence or jump of the specific heat. This is a general result in phase transitions in 
strongly disordered systems and also applies in the present model. An important feature in figure 2 (also observed in 
other studiea!3) is the presence of a maximum in the specific heat at a temperature (~ 3.2) much higher than T c (as 
we will see later, T c ~ 2.62). At T c the specific heat is continuous so we expect the specific heat exponent a to be 
negative. 
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V. THE ORDER PARAMETER 



It was shown long time ago by Edwards and Anderson that the appropriate order parameter for spin-glasses is a 
measure of the temporal freezing of the local variables (in our case, the spins) of the system. In the general framework 
of spin-glass theory the order parameter is computed through the introduction of replicas in the system. This is the 
natural way of introducing the notion of a distance between two configurations in phase space. We take two replicas 
of the same system {<Ji,Ti} (i.e. two identical Hamiltonians in eq.(Q) with identical realization of the couplings J's). 
Then we define the global overlap Q between the two replicas NQ — J3j=i G i T i an( l evaluate its probability distribution 
Pj{q) averaged over the Gibbs measure ((..)) 

Pj(q) = (S(q - Q)> ■ (9) 

In our specific model the overlap is NQ — Yli=i( a i T l + a 2 T 2) where a\, a\ and r|, r\ occupy the site i in the two 
different replicas respectively. Pj{q) gives the probability that two equilibrium configurations <jj, 73 have an overlap q. 
According to the mean-field scenario (the validity of which we would like to check for the present model) ergodicity 
breaks at low temperatures and the phase space splits up into a large number of single ergodic components or states. 
The barriers separating these components diverge with the size of the system suggesting that a symmetry is broken. 
This symmetry is generally referred to as replica symmetry and it is the symmetry under the group of permutations of 
a finite number of replicas. Somehow, this symmetry is artificial (actually, it emerges from the use of the replica trick, 
a general method to deal with the averaging of the logarithm of the partition function in disordered systems). But 
its physical meaning is quite appealing. Different equilibrium configurations <7i,Ti can take different values according 
to the basins of attractions (corresponding to different ergodic components) to which they belong. The function 
P.j(q) is highly non trivial and this is a signature that different configurations quite far one from each other in the 
phase space contribute with a finite weight to the equilibrium partition function. Consequently, different states have 
always the same free energy, internal energy and entropy per site but they differ only in the structure of their typical 
configurations. This is signalled by a non-trivial Pj{q) distribution (i.e. with several peaks at different values of q). 

Another general consequence of the splitting up of the phase space into different states is the possible existence of 
chaotic effects in the equilibrium phaseESa. A small perturbation in the Hamiltonian can change the shape of the 
states as well as reshuffle their Boltzmann weights. After perturbing the Hamiltonian the new configurations can be 
very different from the initial ones. One consequence of this effect is the existence of non self-averaging quantities. 
In particular, if we change completely the microscopic realization of the disorder in the original Hamiltonian (for 
instance, by changing the couplings J in eqs. (^,^)Ell) the new equilibrium states differ completely from the previous 
ones. In this case we do not add energy to the system (the new and the old states have always the same energy per 
site) but the reshuffling of the Boltzmann weights of the different states is enough to change completely the form of 
the Pj(q). Then, the Pj{q) is strongly non-self averaging, a result which has been proved in mean-field theory and 
which we would like to check also in short-range systems. 

The purpose of this section is to show how the study of the P(q) = Pj(q) averaged over the disorder can yield 
evidence for a phase transition in the present model. In the next section we will show that the non self-averaging 
character of the Pj(q) can be used as an independent check for the transition. 

A good way to characterize the P(q) is through its moments. In particular, the first moment g, the second cumulant 
which directly yields the spin-glass susceptibility XSG, the sk ewness Y and the Binder parameter Z. More precisely, 
if we define the average [f(q)] = J dqf(q)P(q) (where P(q) = Pj(q)) then we have : 

XSG 

Y 

* = i(3 



V[(q-[q}) 2 } 

l(q-[q}) 3 } 
l(q-[q}) 2 }* 
[(<?-M) 4 K 

l(q-[q]) 2 } 2> 



(10) 
(11) 



(12) 
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FIG. 3. First moment [q] versus temperature (L = 3,4,5,6 correspond to long dashed, short dashed, dot and solid lines). 
Error bars are shown for L = 6 

In figures 3,4,6,7 we show [q], xsd Y and Z as a function of temperature for different sizes. In figure 3 we show 
the first moment as a function of temperature. In the presence of time reversal symmetry, [q] vanishes (as well as the 
skewness Y) but not in the present case. Note that the curve of [q] as a function of temperature is smooth without 
any sign of a jump or discontinuity. In mean-field spin glass transitions (continuous or discontinuous) this result is 
expected because the jump in [q] is proportional to (1 — m)qEA where qEA is the Edwards- Anderson parameter (the 
maximum value of q such-jthat P{q) ^ 0) and m is the replica symmetry breaking parameter (the size of the blocks in 
the breaking ansatz, seeE3 for general introductory textbooks). In continuous replica symmetry breaking transitions 
qEA vanishes at T c but in first-order replica symmetry breaking transitions qEA is finite at T c and m(T c ) = 1. In both 
cases there is no jump in [q] . 
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FIG. 4. xsg versus temperature. From bottom to top, L — 3, 4, 5, 6. 

The results for the spin-glass susceptibility are more interesting. Experimentally, spin glasses show a divergence of 
the non-linear susceptibility Xni defined through the expansion 
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M(H)= X0 H + XniH 3 + O(H 5 ) . (13) 

It can be generally shownE] that the non-linear susceptibility is related to the spin-glass susceptibility defined in 
cq.([l(]). Although the linear coefficient \o hi the expansion (|l^) does not show any indication of T c the behavior of 
the non- linear term \ni is singular at T c . In mean-field spin glasses with continuous transition xsg shows a power 
law divergence at T c . Contrarily, in mean-field spin glasses with discontinuous RSB xsg shows a finite jump at T c . 
This last feature has been claimed to be,the explanation for the violation found in experiments at the glass transition 
for one of the two Ehrenfest relationsca Below T c , xsg is infinite in both cases. This result is related to the non 
trivial character of the P(q) which has contributions from different values of q. In figure 4 we show XSG for different 
sizes as a function of T. Indeed our results in figure 4 show an algebraic divergence of the spin-glass susceptibility 
Xsg and a least-squares fit of the data in the high-temperature region (where finite size effects are negligible) yields 
XSG ~ (T - T c ) - T with T c ~ 2.63 and 7 ~ 1.0. A finite size scaling plot of the data for xsg/L 2 ^ ~ (T - T c )L l l v is 
shown in figure 5 with v ~2/3, rj ~ 1/2. This is in agreement with the exponents relation 7 = (2 — t])v and with the 
previously estimated value of 7. Using the hyperscaling relation a = 2 — Dv we get a ~ —2/3 in agreement with the 
absence of any singularity or jump in the specific heat (see figure 2). An independent measure of v will be obtained 
in the following section. Anyway, from this first data, we may conclude that the divergence of xsg is related to the 
divergence of a correlation length at T c . 
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FIG. 5. Finite size scaling of the spin-glass susceptibility. Our data is compatible with r\ ~ 1/2, v — 2/3. 
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FIG. 6. Skewness versus temperature. From top to bottom, L = 3, 4, 5, 6. 
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FIG. 7. Binder parameter versus temperature. From bottom to top in the high temperature phase, L = 3, 4, 5, 6. 

Figures 4 and 5 are the first evidence for a phase transition in the model. Figures 6 and 7 show the skewness and the 
Binder parameter as a function of temperature. Because of the adimensional character of these quantities one expects 
that they should be universal in the critical point and related to the amplitudes of the renormalization group flow 
equations. In the large volume asymptotic regime the value of the skewness Y and the Binder parameter Z should be 
volume independent at T c . Consequently T c should manifest as a common crossing point of the curves corresponding 
to different system sizes. A common crossing point can be hardly identified in figures 6 and 7. Nevertheless, the fact 
that Y and Z do not vanish at low temperatures is a sign for a non trivial low temperature phase. Actually these two 
figures yield few information about the transition and it is hard to guess what is the character of thfi-phase transition. 
Let us note that strikingly similar results were obtained for the Ising spin glass in a magnetic fielco. 









0.2 0.4 0.6 

T=2. ,2.353,2.66,3.077,4. for L=6 
FIG. 8. P{q) for L = 6 at different temperatures (from left to right, T = 4, 3.077, 2.66, 2.353, 2.0). 

Gontrarily to what is found in models with time-reversal symmetry the P(q) is not symmetric around q = 0. This 
explains why [q] and the skewness (Y) are non zero. Moreover, the fact that Y < is related to the asymmetric 
character of the P(q). This is clearly shown in figures 8,9 and 10. Figure 8 shows the P(q) for the largest simulated 
size L = 6 at different temperatures T = 4.0, 3.077, 2.66, 2.353, 2.0. Note that above T c the P(q) is Gaussian and 
develops a non-trivial shape at low temperatures with a non vanishing tail which extends down to values of q close to 
zero (the presence of this tail is clearly appreciated plotting the vertical axis in a logarithmic scale). 

pExcept for the fact that the P(q) is not symmetric around q — this behavior is reminiscent of what is observed 
inllj for the p = 4 case. Figures 9 and 10 show the P(q) for different sizes at T = 2.66 and at T = 2.0 respectively. 
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T = 2.66 

FIG. 9. P(q) at T ~ Tc for L = 3, 4, 5, 6. For the sake of clarity, we show the error bars only for L = 6 and only for some 
values of q. 
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-0.5 0.5 1 

T = 2. 

FIG. 10. P(q) at T = 2.0 for L — 3, 4, 5, 6. Errors are shown for L — 6, for some vales of q. 

VI. A NEW EVIDENCE FOR THE TRANSITION 

A new evidence for the existence of a phase transition can be obtained from the study of the sample to sample 
fluctuations of the order parameter function Pj(q). Our main observation originates from recent results obtained by F. 
Guerrao. Guerra has shown that sample to sample fluctuations of the cumulants of the order parameter distribution 
Pj{q) are Gaussian distributed in the thermodynamic limit. Let us now define a sample dependent (i.e. J dependent) 
susceptibility through 

X S G = V({q 2 ) (qf) . (14) 



Note that Xsg * s different from xsg defined in eq.(|iC)|). It can be proveda that the following relationship is fulfilled 
in spin glasses below T c 



n = ix '^ )2 ~ ^ o^_ = i (15) 

vm q -(q)r)-(xi G ) 3 

where, as before, (.) means average over the quenched disorder. 

The interest of defining the parameter G is that it vanishes above the transition temperature in the disordered 
phase where sample to sample fluctuations of the Pj(q) disappear in the V — > oo limit. Similar information to that 
obtained from (|l|) can be also gathered from the sample to sample fluctuations of Xsg > 



A = ^sg) 2 _J^sg) (16) 
(x J sg) 2 

In principle, eq. ( |l6|) yields also non trivial behavior in the low temperature phase even though (in contrast to G) it 
does not necessarily converge (in the thermodynamic limit) to a temperature independent value. Both parameters (A 
and G) are good indicators of the transition although only A gives a precise answer to the question whether the order 
parameter is self- averaging or not. The reasoncj is that G may be finite even when the numerator and denominator jti 
eq.(|l^) vanish. Actually A is the numerator of G so it gives a precise information whether self-averaging is satisfiedE3. 

Note that G is a parameter which plays the same role as the usual Binder parameter g in ferromagnets and is given 
(in the V — > oo limit) by G(T) = (1/3)(1 — Ojj(T — T c )) where Qh is the Heaviside theta function. In RSB transitions 
( |l5| ) goes to zero (as the size V increases) as 1/V for T > T c but converges to a finite value for T < T c . We expect 
the critical temperature (where RS breaks) to be signaled by the crossing of the different curves corresponding to 
different lattice sizes. 
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FIG. 11. Parameter G for L — 3,4,5,6 (open circles, triangles, diamonds and crosses respectively). 

Our results for G and A are shown in figures 11 and 12. Both figures show essentially the same result, i.e. the curves 
for G and A for different sizes display a common crossing point located approximately at T c ~ 2.63 in agreement 
with the result derived in the previous section from the divergence of the spin-glass susceptibility. Assuming for the 
parameter G the following scaling behavior G{T) = G{L/£) with £ - (T-T c )- U then {dG/dT) T=Tc ~ L 1 '" . In figure 
13 we show the scaling behavior for A and G. The scaling plots for A and G yield a more precise fit to the critical 
exponent v because it does involve only one free parameter (T c was obtained looking at the crossing of the different 
curves). A good estimate for v yields v ~ 2/3 for both A and G but precision is not good enough to exclude a slightly 
smaller value (such as v — 1/2). The value 2/3 is in good agreement with the one obtained from the divergence of 
the spin-glass susceptibility. 
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FIG. 12. Parameter A for L — 3, 4, 5, 6 (open circles, rhombi, triangles and crosses respectively) 
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FIG. 13. Parameters A, G for different sizes L = 3, 4, 5, 6 versus (T - 2.62)L 3/2 
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Figures 11 and 12 clearly suggest that replica symmetry breaking takes place below T c . This supports the result 
that the breaking of ergodicity is intimately related to the non self-averagi»g character of the order parameter. This 
result is in contradiction with heuristic arguments by Newman and SteincZI who have suggested that self-averaging 
should be automatically restored in short-range systems due to the translational invariance symmetry of the lattice. 

VII. THE DYNAMICAL EXPONENT Z 

Now that we have corroborated the existence of a thermodynamic phase transition in the model we would like to 
learn more about its nature. In particular we would like to clarify whether the relaxation time (which is the analogous 
of the shear viscosity of real glasses) shows an anomalous behavior in the vicinity of the glass region. It is well known 
that an activated behavior in the relaxation time is one of the main characteristics in real glasses. On the other hand, 
a power law divergence of the relaxation time is a signature for a second order phase transition where a massive mode 
vanishes. 

We have computed the equilibrium time correlation function C(t) at several temperatures above the estimated T c . 
C(t) is defined through 

i V 

c ^ = v E^ )^) + 4(oH(t)) ■ (17) 

i=l 

Compared to similar studies undertaken in the p = 4 caseBIll, the analysis in the present case turns out to be 
more difficult. This is due to the fact that in the absence of time reversal symmetry in the Hamiltonian, C(t) does 
not decay to zero and there is one more unknown parameter (C(oo)). To have a reasonable estimate of z we did two 
types of measures. On the one hand we measured C{t) and fitted it to a stretched exponential form 

C(t) = q+(2-q)cM-(t/r) t3 ) , (18) 

with three free parameters q (C(oo)), t (the relaxation time) and j3 (the stretching exponent). Note that C(t) is 
normalized such that C(0) = 2. Figure 14 shows some of the fits which turn out to be quite good. In this way we 
got some estimates for r which unfortunately are not very precise to determine the value of z. To be more precise 
one should include a power law term t~ a multiplying the exponential terrn-in the fitting function (|l^) as was done by 
Ogielsky in the study of the three dimensional Edwards- Anderson modeE3. Unfortunately including a term of this 
type in ( [i"8| ) introduces too many free parameters into C{t) making fits poorly predictive. Nevertheless, from figure 
14, we may conclude that relaxation turns out to be very slow close to T c 

Our values estimated for the relaxation time exclude any activated behavior. This excludes the existence of a 
viscosity anomaly as well as the existence of two step-celaxation processes in this model. The same conclusion was 
reached in the p = 4, D = 3 case by studying the C(t)t3. _ 

An estimate for z can be obtained by studying the off-equilibrium decay of the order parameter^ or the internal 
energyEl This last case has been applied also to the study of structural glass modelsEJ as well as in the p = 4 caseO. 
In this case, one studies the decay of the internal energy starting from a random initial configuration at T = T c and 
using a fit to a power law behavior of the following form 

E(t) = E{oo) + Ar x . (19) 

Note that this is an off-equilibrium measure which is expected to yield the equilibrium dynamical exponent. Under 
the assumption that hyperscaling is valid, and using simple scaling relations one obtains the exponents relation 
A = (d — \ jv)j z. In figure 15 we show the decay of the internal energy at the estimated T c . A good fit is obtained 
with A = 0.55±0.1 which yields z ~ 4.5± 1. This value for z is very similar to the one found in the Edwards- Anderson 
model in four dimensions. Still, in this model the divergence of the relaxation time is not fast (zv ~ 3) if compared 
with that of the the M = 3,p = 4 modelE3E3 (z ~ 7, zv ~ 6 ± 1). The value of zv being not very large in our model 
(at least, compared to those generally found in spin-glass models in three dimensions) explains why we succeeded in 
getting very clean results through finite-size scaling for the existence of a phase transition. Thermalization was easier 
to achieve. 
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FIG. 15. Energy versus t' - 55 at T = 2.62 for one sample with L = 20 
VIII. DISCUSSION 

In this paper we have investigated the critical behavior of a three spin model in finite dimensions. The motivation was 
to understand how short-range effects modify the AGM mean-field scenario for the spin-glass transition. Moreover the 
Hamiltonian for the present model [M — 2,p = 3) has no time-reversal symmetry. In the presence of a thermodynamic 
transition this lack of time reversal symmetry has far reaching consequences on the nature of the transition. 

We have established (through finite size scaling methods) the existence of a phase transition without latent heat 
and with an algebraic divergence of the spin-glass susceptibility. On the other hand, we find indications that the 
relaxation time diverges according to an algebraic power law as in ordinary continuous phase transitions excluding 
the presence of an activated relaxation time. 

Consequently, we are lead to the conclusion that the first order character of the transition present in the mean-field 
limit (M — > oo) is lost in finite dimensions. In particular, our results v ~ 2/3, r\ ~ 1/2, 7 ~ 1, z ~ 4.5 yield reasonable 
fits to all the data. Moreover, these exponents yield a negative value for the specific heat exponent in agreement with 
the fact that there is not jump or divergence in the specific heat at T c . Note that a vanishes in mean-field so the 
main effect of finite-dimensional corrections is to decrease the value of a. A negative value for a was-apparently also 
obtained for p = 4 in three dimensions for M = 3, 4. Although large size simulations in p = 4, D = 3E3 yield a smaller 
value of v (compatible with v = and hence a = 0) we must exclude this possibility from the absence of any jump 
in the specific heat at T c . I— . 

Quite long ago Gross, Kanter and Sompolinsky, from the solution of the mean-field Potts glasala, suggested the 
possibility that v — 2/D could be valid in finite dimensions similarly as happens for pure systems (with the corre- 
sponding relation v — 1/D). In particular, the explanation for the rounding of the phase transition due to finite size 
effects would be very similar to the explanation valid in first order transitions in pure systemsHcj but now modified 
to account for the presence of randomness. If the transition in finite dimensions were first order then we would expect 
the validity of the relation v = 2/D as well as the absence of upper critical dimension. Our numerical results tend to 
discard this possibility. 

Note that the model we are considering here has no time-reversal symmetry. Consequently, any thermodynamic 
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transition cannot be associated to the breaking of an original symmetry of the Hamiltonian. In this respect, the 
transition we are facing reminds a lot of the mean-field transition of spin-glasses in a magnetic field. We believe that 
the type of transition presented here is one of the most clear examples of second order phase transition in strongly 
disordered systems where replica symmetry breaks. Figures 11, 12 and 13 offer good evidence for this result. 

Finally we would like to comment on the behavior of the entropy of the model as a function of the temperature. 
One prominent prediction in the AGM scenario is the collapse of the configurational entropy at T c . Obviously in the 
present model we do not expect that the configurational entropy vanishes at T c since the transition is continuous. In 
figure 16 we plot the total entropy (which is the sum of the configurational and its intrastate part) as a function of (3 
obtained numerically by integrating the internal energy between (3 = and [3. As (3 increases (data is shown between 
(3 = and /3 = 0.32) the entropy becomes steadily linear and seems to extrapolate to zero at a finite value. A linear 
fit to the last set of data yields S((3) ~ 5.78(0.364 — (3) which vanishes at T c ~ 2.75, a result strikingly close to the 
previously estimated value of T c through finite-size scaling methods. This value of T c is quite stable against data 
inclusion or exclusion in the fit. Actually the same linear behavior is observed for the entropy when plotted versus 
the temperature. A linear fit yields in that case, S(T) ~ 0.536(T — 2.64) a result still very close to previous estimates 
of T c . Consequently the total entropy (and also the configurational entropy) seems indeed to vanish at a temperature 
slightly above T c . This is non sense because above T c the entropy must always be finite. We conclude that the entropy 
must stop decreasing at temperatures close to T c and depart from the linear behavior again. Actually, this is what 
we expect from the presence of the maximum of the specific heat (figure 2). Using the relation C = we safely 
predict a breakdown of the linear behavior at T ~ 3.2 {(3 ~ 0.31) where the entropy should start to form a plateau. 

To conclude, we have found that short-range effects in some class of models (with non-translational invariant 
disorder) make the transition to become second order. The type of transition should be the same to the one expected 
for the Edwards- Anderson model in a magnetic field. The main difference is that, in the last case, the parameter space 
contains the temperature and the field while in the former case the only parameter which controls the phase transition 
is the temperature. This implies strong crossover effects in the critical region for the Edwards- Anderson model in a 
field due to the proximity of the zero-field fixed point. Such crossover effects are not present in the present model 
making the determination of the critical behavior much simpler. It would be very interesting to extend the research of 
strongly disordered spin models without time-reversal symmetry to other cases, such as models with translationally 
spatial invariant disorder, to understand under which conditions the first order character of the mean-field transition 
survives in finite dimensions. 
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FIG. 16. Entropy of the model versus temperature. The continuous line is a linear extrapolation of the last set of points. 
Data were obtained simulating a sample with L = 10. The arrow corresponds to the estimated value of T c . 
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